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We decompose U(l) gauge fields into a monopole and photon part across the phase 
transition from the confinement to the Coulomb phase. We analyze the leading 
Lyapunov exponents of such gauge field configurations on the lattice which are 
initialized by quantum Monte Carlo simulations. It turns out that there is a strong 
relation between the sizes of the monopole density and the Lyapunov exponent. 



1 Motivation 

The main question in QCD is the mechanism of confining quarks into hadrons. 
One seeks for special classes of gauge fields which are responsible for confine- 
ment. In the framework of the dual superconductor, color magnetic monopolcs 
arise as the appropriate candidates. Their investigation necessitates a projec- 
tion of the non-Abelian gauge theory on the compact U(l) theory. In lattice 
computations it was demonstrated that monopoles account for more than 90 
percent of the string tension in the confinement and they drop toward zero 
across the phase transition. 

The study of chaotic dynamics of classical field configurations in field the- 
ory finds its motivation in phcnomcnological applications as well as for the 
understanding of basic principles. The role of chaotic field dynamics for the 
confinement of quarks is a longstanding question. Here, we analyze the leading 
Lyapunov exponents of compact U(l) configurations on the lattice. The real- 
time evolution of the classical field equations was initialized from Euclidean 
equilibrium configurations created by quantum Monte Carlo simulations. This 
way we expect to see a relationship between the appearance of monopoles and 
the strength of chaotic behavior in lattice simulations. 

2 Monopoles in compact quantum field theories 

We have investigated 4d U(l) gauge theory described by the action 

S{*7;}=$>-cos0 p ), (1) 
p 

with Ui = U Xlfl = exp(i6 x<fl ) and 9 p = 6 Xjll + 9 x+ jx. v - 6 X+Dtlx - 6 X , V (v ^ fi) . 
At critical coupling (3 C « 1.01 U(l) gauge theory undergoes a phase transition 
between a confinement phase with mass gap and monopole excitations for 
(3 < [3 C and the Coulomb phase with a massless photon for (3 > (3 C , see Fig. 1. 
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Figure 1: Average of the absolute Polyakov loop squared (|VK| 2 ) as a function of j3 on a 12 4 
lattice. 

We are interested in the relationship between monopoles of U(l) gauge 
theory and classical chaos across this phase transition. Following Refs. [I], |[ 
|3|, |], we have factorized our gauge configurations into monopole and photon 
fields. The U(l) plaquette angles 9 X ^ V are decomposed into the "physical" 
electromagnetic flux through the plaquette 9 X ^ V and a number m x ^ v of Dirac 
strings passing through the plaquette 

9 x .jzv — 9 X ^ v -\- llT 7Yl x ^y , (2) 

where 6 X ^ V E (— tt, +tt] and m x ^ v ^ is called a Dirac plaquette. Monopole 
and photon fields are then defined in the following way 

0£° n = -27r^G x , x ,d> x>;i (3) 

x' 

. (4) 

x' 

Here d' acts on x' , the quantities m x ^ v and 9 X ^ V are defined in Eq. (|J) and 
G x . x i is the lattice Coulomb propagator. One can show that 9 X ^ = 9™™+9v h ° t 
is up to a gauge transformation identical with the original 9 Xi{i defined by 
U Xlll = exp(i9 x ^). 

3 Classical chaotic dynamics from quantum Monte Carlo initial 

states 

Chaotic dynamics in general is characterized by the spectrum of Lyapunov 
exponents. These exponents, if they are positive, reflect an exponential diver- 
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gence of initially adjacent configurations. In case of symmetries inherent in the 
Hamiltonian of the system there are corresponding zero values of these expo- 
nents. Finally negative exponents belong to irrelevant directions in the phase 
space: perturbation components in these directions die out exponentially. Pure 
gauge fields on the lattice show a characteristic Lyapunov spectrum consisting 
of one third of each kind of exponents. □ Assuming this general structure of the 
Lyapunov spectrum we investigate presently its magnitude only, namely the 
maximal value of the Lyapunov exponent, L max . 

The general definition of the Lyapunov exponent is based on a distance 
measure d(t) in phase space, 

L := lim lim —In ,, \ . (5) 

t^ood(o)^o t d(Q) 

In case of conservative dynamics the sum of all Lyapunov exponents is zero 
according to Liouville's theorem, Li — 0. We utilize the gauge invariant 
distance measure consisting of the local differences of energy densities between 
two 3d field configurations on the lattice: 

d: =^E P \trUp-tTU' P \. (6) 

Here the symbol stands for the sum over all Np plaquettes, so this distance 
is bound in the interval (0,27V) for the group SU(N). Up and U' P are the 
familiar plaquette variables, constructed from the basic link variables U Xj i, 

U Wti = exp (aA c x ^T c ) , (7) 

located on lattice links pointing from the position x = {x\, Xi,x$) to x + aei. 
The generators of the group are T c = —igr c /2 with r c being the Pauli matrices 
in case of SU(2) and A x i is the vector potential. The elementary plaquette 
variable is constructed for a plaquette with a corner at x and lying in the 

r/-plane as U x ^j = U x sU x+ ijUl + j A ulj. It is related to the magnetic field 
strength B c x k : 

U Xt ij = exp ( , , a r/ /;; , /") . (8) 
The electric field strength E x i is related to the canonically conjugate momen- 
tum P x j = U Xt i via 

(/''',.,/ !,) . (9) 

The Hamiltonian of the lattice gauge field system can be casted into the 
form 



\(P,P) + 1-\(U,V) 



(10) 
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Here the scalar product stands for (A,B) = ^tr(ABt). The staple variable 
V is a sum of triple products of elementary link variables closing a plaquette 
with the chosen link U. This way the Hamiltonian is formally written as a 
sum over link contributions and V plays the role of the classical force acting 
on the link variable U. 

Initial conditions chosen randomly with a given-average magnetic energy 
per plaquette have been investigated in past years. □ In the present study we 
prepare the initial field configurations from a standard four dimensional Eu- 
clidean Monte Carlo program on a 12 3 x 4 lattice varying the inverse gauge 
coupling /3.U We relate such four dimensional Euclidean lattice field configura- 
tions to Minkowskian momenta and fields for the three dimensional Hamilto- 
nian simulation by identifying a fixed time slice of the four dimensional lattice. 

4 Chaos and confinement 

We start the presentation of our results with a characteristic example of the 
time evolution of the distance between initially adjacent configurations. An 
initial state prepared by a standard four dimensional Monte Carlo simulation 
is evolved according to the classical Hamiltonian dynamics in real time. After- 
wards this initial state is rotated locally by group elements which are chosen 
randomly near to the unity. The time evolution of this slightly rotated config- 
uration is then pursued and finally the distance between these two evolutions 
is calculated at the corresponding times. A typical exponential rise of this 
distance followed by a saturation can be inspected in Fig. 2 from an example 
of U(l) gauge theory in the confinement phase and in the Coulomb phase. 
While the saturation is an artifact of the compact distance measure of the lat- 
tice, the exponential rise (the linear rise of the logarithm) can be used for the 
determination of the leading Lyapunov exponent. The left plot exhibits that 
in the confinement phase the original field and its monopole part have similar 
Lyapunov exponents while the photon part has a smaller L max . The right plot 
in the Coulomb phase suggests that all slopes and consequently the Lyapunov 
exponents of all fields decrease substantially. 

We now turn to a comparison with the monopole density in U(l) quantum 
field configurations, 

Qmon = 4y~ ^ / \ m x.fJ.!^\ ■ (11) 

The left plot of Fig. 3 exhibits for a statistics of 100 independent configurations 
that Qmon decreases sharply between the strong and the weak coupling regime. 
It can be seen that the photon fields carry only a few monopoles. 
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Figure 2: Exponentially diverging distance of initially adjacent U(l) field configurations on 
a 12 3 lattice prepared at f3 = 0.9 in the confinement phase (left) and at /3 = 1.1 in the 
Coulomb phase (right). 




Figure 3: Monopole density (left) and Lyapunov exponent (right) of the decomposed U(l) 
fields as a function of coupling. 



The main result of the present study is the dependence of the leading 
Lyapunov exponent L max on the inverse coupling strength /3, displayed in the 
right plot of Fig. 3. As expected the strong coupling phase, where confinement 
of static sources has been established many years ago by proving the area 
law behavior for large Wilson loops, is more chaotic. The transition reflects 
the critical coupling to the Coulomb phase. The right plot of Fig. 3 shows 
that the monopole fields carry Lyapunov exponents of nearly the same size 
as the full U(l) fields. The photon fields yield a non-vanishing value in the 
confinement ascending toward /3 = for randomized fields which indicates that 
the decomposition ([[]) works perfectly for ideal configurations only. 
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5 Conclusion 



We investigated the classical chaotic dynamics of U(l) lattice gauge field con- 
figurations prepared by quantum Monte Carlo simulation. The fields were 
decomposed into a photon and monopole part. The maximal Lyapunov ex- 
ponent shows a pronounced transition as a function of the coupling strength 
indicating that on a finite lattice configurations in the strong coupling phase 
are substantially more chaotic than in the weak coupling regime. The compu- 
tations give evidence that the Lyapunov exponents in the original U(l) field 
and in its monopole part are very similar. The situation for the monopole 
density is analogous and serves as a consistency check of the decomposition. 
We conclude that classical chaos in field configurations and the existence of 
monopoles are intrinsically connected to the confinement of a theory. 
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